Sensitivity of neutron to proton ratio toward the high density behavior of symmetry 

energy in heavy-ion collisions 
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The symmetry energy at sub- and supra-saturation densities has great importance for understand- 
ing the exact nature of asymmetric nuclear matter as well as neutron stars, but it is poorly known, 
especially at supra-saturation densities. We will demonstrate here whether or not the neutron-to- 
proton ratios from different kinds of fragments can determine the supra-saturation behavior of the 
symmetry energy. For this purpose, a series of Sn isotopes were simulated at different incident 
energies using the isospin quantum molecular dynamics (IQMD) model with either a soft or a stiff 
symmetry energy. It is found that the single neutron-to-proton ratio from free nucleons as well as 
Light Charged particles (LCPs) is sensitive to the symmetry energy, incident energy, and isospin 
asymmetry of the system. However, with the double neutron-to-proton ratio, this is true only for 
the free nucleons. It is possible to study the high-density behavior of symmetry energy by using the 
neutron-to-proton ratio from free nucleons. 

PACS numbers: 21.65.Ef, 21.65. Cd, 25.70.Pq, 25.70.-z 
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I. INTRODUCTION 

From the Bethe-Weizsacker mass formula, it is well un- 
derstood that the symmetry energy from bulk matter is 
the difference between the energy of pure neutron matter 
and pure symmetric matter. Mathematically, it can we 
written as 



E Sym {p,S) = E(p,S = l)-E( P ,S = 0), 



(1) 



where 6 = — + p p p and p = p n + p p . p n , p p , and p are the 
neutron, proton, and nuclear matter densities, respec- 
tively. The symmetry energy has great importance in 
the dense matter existing in the neutron stars, but only 
indirect information can be extracted from astrophysi- 
cal observations It is also important in the quark 
gluon plasma (QGP) and hadron gas (HG) phase 0. 
The QGP and HG phases existed in the early stage of 
the evolution of Universe (about 15 billion years ago) 
and are inaccessible nowadays. It is difficult to recre- 
ate these conditions, although numerous experiments are 
occurring at the Relativistic Heavy Ion Collider (RHIC) 
and the Large Hadron Collider (LHC) Q. Heavy- ion 
reactions, during which matter goes through compres- 
sion and expansion, are considered to be the true testing 
ground for the hot and dense matter phases. The nuclear 
equation of state (NEOS) and the density dependence of 
the symmetry energy can be probed by some observablcs 
in intermediate-energy heavy- ion collisions (HICs). The 
softness of the NEOS has been well described in the lit- 
erature in the last couple of decades 4, 5]. However, the 
density dependence of the symmetry energy, from the 



Coulomb barrier to the deconfincment of nuclear matter, 
is a hot topic in the present era. At sub-saturation den- 
sities, the density dependence of the symmetry energy is 
studied by observables such as the neutron-to-proton ra- 
tio, isotopic and isobaric scaling, isospin diffusion, isospin 
fractionation and/or distillation, and isospin migration. 

Recently, the MSU group Q claimed the softness 
of the symmetry energy at sub-saturation densities by us- 
ing the double neutron-to-proton ratio and isospin diffu- 
sion from two isotopic systems, 112 Sn + 112 Sn and 124 Sn 
+ 124 Sn at E = 50 MeV/nucleon. In another study, again 
soft symmetry energy was claimed by using the isospin 
diffusion for the same set of reactions, but at E = 35 
MeV/nucleon Q. In a recent study, soft symmetry en- 
ergy is also favored for the same set of reactions at E 
= 50 MeV/nucleon by using the neutron-to-proton ra- 
tio @. In all the studies, the problem of sub-saturation 
density dependence of the symmetry energy seems to be 
addressed to some extent; however, the uncertainties are 
still large enough to justify the large amount of work that 
is under way in many laboratories all over the world. 

In contrast, the present status of supra-saturation den- 
sity dependence of the symmetry energy is quite uncer- 
tain and interesting. The high-density behavior of the 
symmetry energy in the literature is studied by using 
two important parameters: one is the yield ratio param- 
eter and second is the flow parameter. The yield ratio 
parameter has been studied in term of single and double 
ratios of neutrons to protons ficil - fl2j , single and double 
ratios of 7T~ /7r +JlbT [l3l4l8| , the E"/E+ ratiojli], the 
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K~ /K + ratio [l5|, and isospin fractionation 12j, while, 
the flow parameter has been studied in terms of relative 
and differential flows (single and double ratios) of neu- 
trons to protons or 3 H to 3 He 0, EH, and in terms of 
the ratio [20[ or difference [2l[ of neutron-to-proton el- 
liptic flow. Before using the 3 H and 3 He particle yield 
and flow ratios for the density dependence of the sym- 
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metry energy at high incident energies, one must check 
the production of these particles in the supra-saturation 
density region, which is obtained during the highly com- 
pressed stage only. However, the production of neutrons 
and protons occurs in large amounts and can explain the 
high density dependence of symmetry energy with great 
accuracy, favorable results with neutron and proton el- 
liptic flow at E = 400 MeV/nucleon were also observed 
in 2011. In one of the studies, the softness of the symme- 
try energy with 7^ = 0.9 is predicted by comparing the 
FOPI collaboration data with the neutron-to-proton el- 
liptic flow ratio [2(|. In the same year, Cozma et al. [2l[ 
predicted the softness of symmetry energy with x = 2 
by comparing the FOPI data with the neutron-to-proton 
elliptic flow difference. Even then uncertainty lies in the 
results, in terms of determination of symmetry energy: in 
the first study, symmetry energy is momentum indepen- 
dent, while in later one it is from momentum-dependent 
interactions. Moreover, the studies were limited to only 
400 MeV/nucleon. 

Let us examine some interesting features from the ra- 
tio parameters at supra-saturation densities. All the ratio 
parameters show sensitivity to the symmetry energy. In 
the literature, it is also claimed that K°JK + and 
have more sensitivity than tt~/tt + [2 EH- The sensitiv- 
ities of all the parameters is checked in term of trans- 
verse momentum and rapidity distribution dependence 
EH EH EH s wnue pi° n an d kaon ratio studies are 
extended with the isospin asymmetry of the system and 
the incident energy [13l Il5l4l7j |. In recent years, when 
the pion ratio has been compared with the FOPI data by 
using the two well known models IBUU04 and ImlQMD, 
in terms of isospin asymmetry and incident energy, the 
predictions for the symmetry energy are found to be to- 
tally opposite. ImlQMD predicts stiff symmetry energy 
(7; = 2) fl7|. while IBUU04 predicts soft symmetry en- 
ergy (x=i)[i6[. 

In the present era, the ir~ /tt + ratio is supposed to be a 
strong candidate for predicting the high-density behav- 
ior of symmetry energy. Just as for ir~ and ir + , neu- 
trons and protons are also produced in large amounts 
up to 1 GeV/nuclcon. Even around 400 MeV/nucleon, 
the production of neutrons and protons is greater than 
the production of pions. Unfortunately, the neutron- 
to-proton ratio parameter in most studies is restricted 
only with the transverse momentum and kinetic energy 
dependencies [Iol - [l2"l ]. To draw a fruitful conclusion in 
the near future, first, it is very important to study the 
isospin asymmetry and incident energy dependencies of 
the neutron-to-proton ratio, just as in the recent n~ /ir + 
study, and then compare the sensitivity to the symmetry 
energy from both ratios, as the pion ratio results were 
recently compared with the FOPI experimental findings. 
Second, one has to avoid choosing randomly any type of 
fragment to study the supra-saturation density depen- 
dence of symmetry energy. For this, it is important to 
check whether or not a particular type of fragment is 
formed in the region p > po, which is simple when one 



addresses the sub-saturation density dependence of sym- 
metry energy. Finally, with increasing incident energy, 
the time evolution of the density has a different trend at 
two extremes: one at the time of maximum compression 
and the second at the freeze-out time (t = 200 fm/c). 
It also becomes interesting to see the different stiffnesses 
of the symmetry energy dependence of the neutron-to- 
proton ratio for incident energies at the time of maximum 
compression and at freeze-out time. 

In the concluding remarks of this paper, we have tried 
to address the following goals: 

• To check the sensitivities of different kind of frag- 
ments to the high-density behavior of symmetry 
energy. 

• To check the behavior of the neutron-to-proton ra- 
tio at the time of maximum compression and at 
freeze-out time. 

• To check the isospin asymmetry and incident en- 
ergy dependences of single and double neutron-to- 
proton ratios from different neutron-rich systems to 
the high-density behavior of symmetry energy, and 
then compare the sensitivity of symmetry energy 
from the neutron-to-proton ratio with that from 
the pion ratio. This study is similar to recent stud- 
ies using the pion ratio [l~3|, Iill4l7j . 

For the present study, the isospin quantum molecular 
dynamics (IQMD) model is used to generate the phase 
space of nucleons, which is discussed in Sec. II. The 
results are discussed in sec. Ill, and we summarize the 
results in sec. IV . 



II. METHODOLOGY: ISOSPIN QUANTUM 
MOLECULAR DYNAMICS MODEL 

In the IQMD model [IH [23| , nucleons are represented 
by wave packets, just as in the QMD model of Aiche- 
lin dj]. These wave packets of the target and projectile 
interact via the full Skyrme potential energy, which is 
represented by U and is given as: 

U = U p + Ucoui- (2) 

Here Ucoui is the Coulomb energy and U p originates from 
the density dependence of the nucleon optical potential, 
and is given as 

The first two of the three parameters of Eq. [3] (a and 
/3) are determined by demanding that, at normal nuclear 
matter densities, the binding energy should be equal to 
16 MeV and the total energy should have a minimum at 
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Pq. The third parameter 7 is usually treated as a free pa- 
rameter. Its value is given in term of the compressibility: 



K = V 



2 



dp 2 \A 



E 



(4) 



The different values of compressibility give rise to soft 
and hard equations of state. The soft equation of state 
is employed in the present study with the parameters 
a = - 356 MeV, @ = 303 MeV, and 7 = 7/6, cor- 
responding to an isoscaler compressibility of K = 200 
MeV. In the third term -Bg°* m is the potential part of 
the symmetry energy, which is adjusted on the basis of 
calculations from the microscopic or phenomenological 
many-body theory, having the form 



ppot 
^Sym 



c s 



(5) 



Here C StP = 35.19 MeV, parameterized on the basis of 
the experimental value of the symmetry energy, is known 
as the symmetry potential energy coefficient. On the 
basis of the 7$ value, symmetry energy is divided into 
two types with 7^ = 0.5 and 7, = 1.5, corresponding 
to soft and stiff symmetry energies, respectively. 

The total symmetry energy per nucleon employed in 
the simulation is the sum of the kinetic and potential 
terms and is given as 



E, 
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25 MeV is known 



as the symmetry kinetic energy coefficient. The kinetic 
symmetry energy originates from the Fermi-Dirac distri- 
bution [24[. 

Finally, we get a density and isospin-single particle po- 
tential in nuclear matter as follows: 



Vr( P ,S) 
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dp 
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(7) 



The 



Here t ^ t', f^- 

' ' dpn P ' OPp P 

potential also depends on the momentum-dependent in- 
teractions, which are optional in the IQMD model. 

Note that the 7 used in the determination of the equa- 
tion of state and ji used in the determination of sym- 
metry energy arc different parameters. The interesting 
feature of symmetry energy is that its value increases 
with decreasing 7, at sub-saturation densities, while the 
opposite is true at supra-saturation densities. In other 
words, soft symmetry energy is more pronounced at sub- 
saturation densities, while stiff symmetry energy is more 
pronounced at supra-saturation densities. 

In the calculations, we use the isospin-dependent in- 
medium cross section in the collision term and the Pauli 



blocking effects as in the QMD model [4|. The cluster 
yields are calculated by means of the coalescence model, 
in which particles with relative momentum smaller than 
PFermi and relative distance smaller than Rq are coa- 
lesced into a one cluster. The value of i?o and PFermi for 
the present work are 3.5 fm and 268 MeV/c, respectively. 



III. RESULTS AND DISCUSSIONS 

The neutron-to-proton ratio is among the first observ- 
ables that was proposed as a possible sensitive probe for 
symmetry energy prediction at sub-saturation densities 
@j 0b however, some studies are also performed using 
the rapidity distribution and transverse momentum de- 
pendencies at supra-saturation densities. In this article, 
the sensitivities of free nuclcons, light charged particles 
(LCPs, having charge number between 1 and 2), and in- 
termediate mass fragments (IMFs, having charge number 
between 3 and ZTotai/6) to the high density behavior of 
symmetry energy are checked, providing the results of 
incident energy and isospin asymmetry dependencies of 
single and double neutron-to-proton ratios with the high- 
density sensitive fragments. 

To perform the study, thousand of events are simulated 
for the isotopes of Sn, namely 112 Sn + 112 Sn, 124 Sn + 
124 Sn and 132 Sn + 132 Sn between incident energies of 50 
and 600 MeV/nucleon at semicentral geometry by using 
the soft and stiff symmetry energies of 7* = 0.5 and 
1.5, respectively. As discussed earlier, a soft equation of 
state with an isospin- dependent nucleon-nucleon (NN) 

cross section of a me d = ^1 — 0-2p~) &free is employed. 
The incident energy and isospin asymmetry dependences 
of single and double neutron-to-proton ratios, just as for 
the n" /n + ratio [HI, 53b are considered as a point of 
importance in the present study. The single ratio is just 
the ratio of neutrons to protons and is represented in 
the study by R(N/Z), while double ratio is the ratio of 
the single ratios of any two isotopes of Sn. In order to 
study the systematics of the isospin effects, the single 
ratio of the isotope with a greater number of neutrons 
is always mentioned in the numerator when the double 
ratio is calculated. Mathematically, the double ratio is 
represented by DR(N/Z) and is given as 



DR{N/Z) = 



R(N/Z) 



neutron rich 



fifjy I %\n>eutron weak 



(8) 



To predict the high-density behavior of symmetry en- 
ergy, the very first point is to understand the time evo- 
lution of the average density at different incident ener- 
gies. With increasing incident energy, the density will be 
expected to be greater than the normal nuclear matter 
density in the most compressed region. As we know, the 
density of the environment surrounding the nucleons of a 
fragment plays a crucial role in determining the physical 
process behind its formation. In Fig. [TJ we display the 
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FIG. 1: (Color online) Time evolution of the average density 
for the 132 Sn + 132 Sn reaction with the soft symmetry en- 
ergy (7^ = 0.5) at semicentral geometry. The different lines 
represent different incident energies ranging from 50 to 600 
MeV/nucleon. 



average density (p/po) reached in a typical reaction as a 
function of time at different incident energies for 132 Sn + 
132 Sn by using the soft symmetry energy ji = 0.5. The 
average nucleon density is calculated as 



(p) = ( 



1 



A T +Ai 



A T +Ap A T +A 

E E 

i=l 3=1 



1 



(2ttL) 3 / 2 
exp[-{fi-f-f/2L]), 



(9) 



with fi and fj being the position coordinates of the i th 
and j th nucleons, respectively. 

As we have expected, with increasing of incident en- 
ergy, the density is found to increase in the compression 
zone. Interestingly, at lower beam energies, the maxi- 
mum density reached is lower and the reaction time is 
longer. With increasing incident energy, the life-time of 
the high-density nuclear matter gets shorter due to in- 
stability. For example, at b = 2 fm the average density 
reaches a maximum and is close to normal nuclear matter 
density at t= 18 and 33 fm/c, respectively, for E = 50 
MeV/nucleon; but for the case of E = 600 MeV/nucleon, 
the respective times are 10 and 20 fm/c. This means 
that the difference between the two times is almost 15 
fm/c at E = 50 MeV/nucleon, while it is only 10 fm/c 
at E = 600 MeV/nucleon. This clearly indicates that the 
matter shows high-density behavior only for a small time 
interval, which decreases with increasing incident energy. 
Since we are interested in the sensitivities of different 
kinds of fragments and their neutron-to-proton ratios, 
only those fragments that lie in the high-density region 
(p > pa) will be sensitive to the high-density behavior of 
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FIG. 2: (Color online) Time evolution of free nucleons (top), 
LCPs (middle), and IMFs (bottom) at semicentral geometry 
for 132 Sn + 132 Sn using the soft symmetry energy (7* = 0.5). 
The different lines have the same meaning as in Fig. [1] The 
vertical line in each panel represents our time limit before 
which the system can be in the supra-saturation region (refer 
to Fig. 1). 



the symmetry energy. 

To check the sensitivities of different kind of fragments, 
in Fig. [2] we display time evolution of free nucleons 
(top), LCPs (middle), and IMFs (bottom) at semicen- 
tral geometry for incident energies ranging from 50 to 600 
MeV/nucleon. The behavior for all kinds of fragments is 
consistent with the results in the literature [25( . The pro- 
duction of free nucleons increases with incident energy, 
and LCP production decreases after 400 MeV/nucleon. 
In Ref. [2j], LCP production is correlated with the nu- 
clear stopping and is also found to have a maximum at 
400 MeV/nucleon. IMF production is found to decrease 
after 100 MeV/nucleon. This is due to the different origin 
of the production of IMFs as compared to free and LCPs. 
For more details about the incident energy dependence 
of IMFs, see Ref. 

Our main task is to check the sensitivities of the frag- 
ments in the high-density region. For this, we apply the 
limit that at least one particle must be produced before 
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FIG. 3: (Color online) Excitation function of (p/po) 7 % which 
is proportional to Es ym , at semicentral geometry. The left 
panel is at the time of maximum compression, while the right 
one is at the freeze-out time. Solid and open circles with 
a dashed line represent the contributions of the soft (•y i =0.5) 
and stiff (74=1.5) symmetry energies, respectively, for the sys- 
tem 132 Sn + 132 Sn. 



the time 20 fm/c, because, in an average, after that time 
the density becomes lower than normal nuclear matter 
density for all the incident energies under consideration. 
The free nucleons are highly sensitive at all the ener- 
gies. This is not true for LCPs and IMFs. LCPs are 
produced in this region only after the incident energy 
reaches 200 MeV/nucleon. In contrast, no IMFs are pro- 
duced in the supra-saturation density region. This means 
that IMFs are not so sensitive to the high-density depen- 
dence of symmetry energy; however, they can be used at 
sub-saturation and saturation densities [7[ . Interestingly, 
the single neutron-to-proton ratio from IMFs is found to 
change with the incident energy (not shown here), but 
this is mainly due to Coulomb interactions. Here we con- 
clude that the neutron-to-proton ratio from free nucleons 
as well as LCPs can act as a probe of the high-density 
behavior of the symmetry energy. 

One more interesting observation is obtained from Fig. 
[TJ With increasing incident energy, the time evolution of 
the density is exactly opposite during the compression 
and expansion stages. That is, in the expansion stage 
the average density is found to decrease with increasing 
incident energy, which was earlier increasing in the com- 
pressed zone. Now, we have two aspects of the basis of 
the time evolution of density: one is the compressed-zone 
time and second is the freeze-out time. Interestingly, if 
the density behavior is opposite at the two times, then it 
would supposedly affect the magnitude of the symmetry 
energy as well as its effect on the nuclear matter during 
the whole time evolution. 

To see the virtual change in the symmetry energy due 
to the change in the density, we display in Fig. [3] the 
incident energy dependence of (p/po) 7 S which is propor- 



tional to the symmetry energy, at the time of the maxi- 
mum compression (left panel) and at the freeze-out time 
(right panel) . At the time of maximum compression, the 
symmetry energy rises with the increasing incident en- 
ergy (increase in density). As the density is more than 
the normal nuclear matter density in this region, the stiff 
symmetry energy is stronger than the soft one. With 
increasing incident energy (increase in density), the stiff 
symmetry energy is changing drastically, while, the soft 
symmetry energy shows little change. This exactly co- 
incided with the ideal picture of density dependence of 
the symmetry energy. On the other hand, if we look 
at the energy dependence of (p/po) 11 at i = 200 fm/c, 
the situation is totally different. The symmetry energy 
is found to decrease with increasing incident energy (de- 
crease in density). Now the density is lower than the 
normal nuclear matter density, so the magnitude of the 
soft symmetry energy is greater than that of the stiff 
symmetry energy. In other words, the supra-saturation 
(sub-saturation) density region is more neutron rich with 
7t = 1.5 (-fi = 0.5). The effect from the sub- and 
supra-saturation density behaviors of symmetry energy 
will compete and contribute in the final obscrvables. Due 
to the different behavior of density at different times, it 
is important to observe the isospin effects at the time of 
maximum compression and at the freeze-out time to un- 
derstand the high-density behavior of symmetry energy. 
For this purpose, in the coming sections, the incident en- 
ergy and isospin asymmetry dependencies of the single 
and double neutron-to-proton ratios from free nucleons 
and LCPs are analyzed. 



A. Incident energy and isospin asymmetry 
dependencies of single neutron to proton ratio 

In order to address the sensitivity of symmetry energy 
at the time of maximum compression and at freeze-out 
time, we display the incident energy and isospin asym- 
metry dependencies of the single ncutron-to-proton ratio 
at different times in Figs. 21 [S] and [5J The left and 
right panels are at the time of maximum compression 
and freeze-out time, respectively. In Fig. [4j the ratios 
from free nucleons and LCPs are displayed in the top 
and bottom panels. The many interesting facts are re- 
vealed in the figure. The incident energy dependencies of 
the ratios are found to be highly sensitive to symmetry 
energy for the two different times. As we know, the rela- 
tive strength of symmetry energy is opposite at sub- and 
supra-saturation densities with 7* = 0.5 and 7$ = 1.5. In 
the range of 50-150 MeV/nucleon, only the low density 
part up to about l.lpo contributes. Therefore, in the 
low-energy region, for free nucleons, we can see the high 
ratio with the soft symmetry energy at both the times 
under consideration. At and above 200 MeV/nucleon, 
a broad range of densities up to 1.8po is involved. Of 
course, at about 200 MeV/nucleon and above, for the 
behavior of the high-density symmetry energy, there is a 
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FIG. 4: (Color online) Excitation function of the single 
neutron-to-proton ratio for free nucleons (top panel) and 
LCPs (bottom panel) at semicentral geometry. The left and 
right panels are the same expect they are at the time of max- 
imum compression and at the freeze-out time, respectively. 
Solid and open circles represent the soft and stiff symmetry 
energies, respectively. Solid and dashed lines corresponds to 
the systems of 112 Sn + 112 Sn and 124 Sn + 124 Sn, respectively. 



combined effect for particles going through both low-and 
high-density region. At higher energies, a higher N/Z 
ratio is observed with the stiff symmetry energy for the 
neutron rich system 132 Sn + 132 Sn at the time of maxi- 
mum compression, which is true with the soft symmetry 
energy at the freeze-out time. The result is similar for 
free nucleons and LCPs. However, LCPs are not as sensi- 
tive and the ratio is even less than the ratio of the system. 
This is due to the excess number of protons involved in 
the production of LCPs compared to free nucleons. These 
protons will lower the ratio for LCPs. 

It is clearly visible that the ratio at both times is almost 
the same with the stiff symmetry energy, but changes 
drastically with the soft symmetry energy. This is due 
to the fact that, at the time of maximum compression, 
the density is in the supra-saturation region and the stiff 
symmetry energy is much higher (see Fig|3]) than the soft 
symmetry energy. Therefore, the stiff symmetry energy 
is able to separate most of the neutrons near the time 
of maximum compression and then accelerate the neu- 
trons toward higher kinetic energy at later times. How- 
ever, the soft symmetry energy is not so high, and the 



separation of the neutrons takes place for a longer time. 
After 50-60 fm/c (see Fig. [1]), the density drops to the 
sub-saturation density region and now the soft symmetry 
energy has a quite high magnitude (see Fig. [3]) compared 
to the stiff one. The soft symmetry energy in this region 
is still separating the neutrons as well as accelerating 
them high kinetic energy. That is why the ratio with the 
soft symmetry energy drastically changes when one goes 
from compression to freeze-out time, but remains almost 
constant with the stiff symmetry energy. 

Mainly, the neutron-to-proton ratio is found to de- 
crease with the incident energy for free nucleons as well 
as for LCPs, just like the ir~ /ir + ratio. The decrease in 
the ratio may be due to two reasons: 

• One reason may be the role of Coulomb interac- 
tions with incident energy. With increasing inci- 
dent energy, chances of break-up of initial corre- 
lations among the nucleons becomes stronger, and 
the production of free nucleons including neutrons 
and protons will increase . However, at very low 
incident energy, the production of neutrons is more 
due to the symmetry energy because of its repul- 
sive (attractive) nature for neutrons (protons). In 
short, due to Coulomb interactions, a shift of pro- 
tons takes place from low to high incident ener- 
gies. The effect of the Coulomb interactions can be 
checked by taking the double ratio, which is dis- 
cussed in Fig. 

• The contribution of pions from secondary-chance 
nucleon-nuclcon collisions might increase with the 
beam energy. If a first-chance nuclcon-nucleon col- 
lision converts a neutron to a proton by producing 
a 7r~, then subsequent collisions of the energetic 
protons can convert them back to neutrons by pro- 
ducing a 7r + . Therefore, at sufficiently high energy, 
the neutrons, which are produced due to symme- 
try energy, are changing into the protons and fur- 
ther producing it's, which will lead to a decease 
in the neutron-to-proton ratio. This can be con- 
firmed by using the double ratio concept. If the 
double ratio is still deceasing with incident energy, 
then it means that, in addition to the Coulomb in- 
teractions, the phenomenon of secondary nucleon- 
nucleon collisions is also very important. 

One more point of interest is that the difference be- 
tween the soft and stiff symmetry energies at freeze-out 
time is found to decrease with incident energy for free 
nucleons, while it increases for LCPs. Of the above two 
reasons, the first one is applicable for free nucleons as 
well as for LCPs. The second one is applicable only for 
free nucleons, as the energy in this study is up to 600 
MeV/nucleon, which is quite sufficient to produce pions. 

To see the effect of the high-density behavior of symme- 
try energy on the isospin asymmetry dependence, we dis- 
play the ratio from free nucleons and LCPs in Figs. [S]and 
[S]at only high energies (200, 400, and 600 MeV/nucleon). 
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Due to the instability of the highly compressed zone, we 
are not able to differentiate between the results of sym- 
metry energy obtained at the time of maximum compres- 
sion; however, we had earlier obtained some important 
conclusions from Fig. |4l where incident energy depen- 
dence was discussed. 

The results from Figs. [5] and |5] at the freeze-out time 
reveal many important points. The isospin asymmetry 
dependence of the ratio from free nucleons is highly sen- 
sitive to the symmetry energy compared to LCPs, i.e., the 
ratio from free nucleons is found to be sharply increas- 
ing with the isospin asymmetry of the system compared 
to LCPs. This is due to the fact that isospin effects on 
the ratio from free nucleons are strongly affected by the 
symmetry energy and weakly affected by Coulomb in- 
teractions, while the opposite is true for the ratio from 
LCPs. As discussed earlier, the ratio is found to decease 
with the incident energy, which is also true here for the 
isospin asymmetry dependence. 

The difference between the soft and stiff symmetry en- 
ergy results comes from the behavior of free nucleon emis- 
sion with the isospin asymmetry of the system, i.e., the 
greater the isospin asymmetry of the system, the greater 
the contribution of neutrons in the ratio due to the sym- 
metry energy. The soft symmetry energy is stronger at 
the freeze-out time, which will lead to an increase in the 
ratio more sharply than the stiff symmetry energy. This 
effect is again weakly observed in the ratio from LCP's. 
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FIG. 5: Isospin asymmetry dependence of the single neutron- 
to-proton ratio for free nucleons at different incident energies. 
The left panel is at the time of maximum compression, while 
the right panel is at the freeze-out time. Solid and open circles 
represent the soft and stiff symmetry energies, respectively. 



B. Incident energy and isospin asymmetry 
dependencies of double neutron to proton ratio 

In order to cancel the Coulomb effects and to see the 
effect of symmetry energy, we show in Fig. the inci- 
dent energy dependence of the double ratio from different 
isotopes of Sn with different combinations, namely, 132 Sn 
+ 132 Sn and 124 Sn + 124 Sn, 124 Sn + 124 Sn and 112 Sn + 
112 Sn, 132 Sn + 132 Sn and 112 Sn + 112 Sn, having differ- 
ences of 8, 12, and 20 neutrons and the same number 
of protons. The upper and lower panels are for free nu- 
cleons and LCPs. The double ratio is found to increase 
with the difference of the number of neutrons in the dif- 
ferent combinations. It is similar to the results obtained 
at sub-saturation densities by different models @,0|- The 
main point to be discussed here is that the double ratio 
is found to decrease with increasing incident energy. As 
we have discussed in Fig. 01 there may be two reasons for 
the decrease of the single ratio with increasing incident 
energy. One was the Coulomb effect, which is canceled 
out here. The second was the pion effect, which is still 
active in the double ratio and becomes more and more 
dominant with increasing incident energy. Due to that 
effect, the double ratio is found to decrease with the inci- 
dent energy. It indicates that the pion production effect 
is very important at high incident energy and is equally 
useful for understanding the high-density behavior of the 
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FIG. 6: Same as in Fig. [5] but for the LCPs. 
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FIG. 7: (Color online) Excitation function of the double 
neutron-to-proton ratio from different isotopes of Sn for free 
nucleons (top) and LCPs (bottom). The vertical line in 
the bottom panel represent the energy limit above which 
DR(N/Z) of LCPs becomes more or less insensitive. Solid 
and open circles represent the soft and stiff symmetry en- 
ergies, respectively. The solid, dashed, and dot-dashed line 
corresponds to double ratios from 
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124 Sn, 124 Sn + 124 Sn to 112 Sn + 112 Sn, and 
112 Sn + 112 Sn, respectively. 
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FIG. 8: Isospin asymmetry dependence of the double neutron- 
to-proton ratio from free nucleons at different incident ener- 
gies. The different symbols have the same meaning as in Fig. 
® 



neutrons to protons from free nucleons. Another possible 
probe is the it~ /it + ratio, which recently was compared 
with the experimental data of the FOPI by the IBUU04 
and ImlQMD calculations 0, [Hi ■ 



incident energy 0, [I?} . 

In contrast, this effect is valid only for the double ra- 
tio from the free nucleons and not from the LCPs. The 
double ratio from LCPs is found to be constant above 
200 MeV/nucleon. This indicates that the effect of the 
symmetry energy for the ratio from LCPs can be ana- 
lyzed only near sub-saturation densities close to l.lpn- 
The decrease in the single ratio for the LCPs was only 
due to the Coulomb interactions at higher incident en- 
ergies, which is canceled out by taking the double ratio; 
the double ratio from the symmetry energy becomes in- 
dependent of the incident energy after 200 MeV/nucleon. 
This type of dependence for the single 7r~/7r + ratio can 
be observed above 1 GeV/nucleon [l3|. The behavior of 
symmetry energy for the double ratio is exactly the same 
as that for the single ratio. This indicates that LCP pro- 
duction is also not a sensitive probe for investigating the 
high-density behavior of the symmetry energy. The only 
possible probe from the fragments is the double ratio of 



In order to strengthen our conclusion, in Fig. [5J we 
display the isospin asymmetry dependence of the dou- 
ble ratio from free nucleons at different incident energies. 
All the curves arc fitted with a power law of the form 
y = ax T , where y is the double ratio from free nucleons 
and x is the double ratio of the systems. The power-law 
exponent r is found to vary drastically with the symme- 
try energy, which is to be discussed later in Fig. |H1 After 
canceling the Coulomb effects, the trend for the double 
ratio is the same as that of the single ratio in Fig. [5J It 
reflects the fact that the isospin effects for free nucleons 
is stronger for more neutron-rich systems and is mainly 
due to the symmetry energy. However, the decrease in 
the isospin effect with the increase of incident energy is 
due to the production of pions at sufficiently high energy. 
The difference in the double ratio obtained with the soft 
and stiff symmetry energies here is also found to increase 
from the neutron-poor to the neutron- rich system, just 
like the single neutron-to-proton ratio in Fi g. [5] as well 
as the single pion ratio in the literature 1|| 17 1. 
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FIG. 9: Incident energy dependence of the power-law expo- 
nent t from Fig. [5] The symbols and lines are the same as in 
Figs. 5 and 8. 

C. Incident energy dependence of Power law 
exponent r 

To see the clear systematics of the incident energy to- 
ward the symmetry energy, we plot the incident energy 
dependence of the power exponent r in Fig. |Hl which 
is extracted from the curves of Fig. [H With increasing 
incident energy, the sensitivity of the symmetry energy 
goes on decreasing toward the double ratio; however, the 
soft symmetry energy is more sensitive in comparison 
with the stiff one. In brief, when one goes from the sub- 
saturation to the supra-saturation density region, the soft 
symmetry still has a crucial role to play compared to the 
stiff one. This is due to the density (Fig. [T]), which 
undergoes a sudden change between the supra- and sub- 
saturation density regions with time at higher incident 
energies. 

Finally, from this study, we confirm that the high- 
density behavior of symmetry energy can be studied by 
using the single and double ratios of neutrons to protons 



from free nucleons. In comparison, the double ratio is 
more accurate for this purpose, due to its greater sensitiv- 
ity to the soft symmetry energy. Meanwhile, the lighter 
and heavier fragments ratio can be considered good can- 
didates at sub-saturation densities, and also have been 
used in the literature many times by different groups 

sa. 

IV. CONCLUSION 

In order to investigate the high-density behavior of 
the symmetry energy, isospin asymmetry and beam en- 
ergy dependencies of neutron-to-proton ratios (single and 
double) from different kinds of fragments are studied by 
using the IQMD model. The single neutron-to-proton 
ratio from free nucleons and LCPs is found to decrease 
(increase) with incident energy (with the isospin asym- 
metry of the system). Stronger isospin effects are ob- 
served with the soft symmetry energy. Similar results 
with the ir~ /tt + ratio are also observed by Li et al. and 
Feng et al., but with opposite behavior for symmetry 
energy. The double neutron-to-proton ratio from free 
nucleons is highly sensitive to the symmetry energy, inci- 
dent energy, and isospin asymmetry of the system. How- 
ever, the sensitivity of the neutron-to-proton double ra- 
tio from LCPs to the nuclear symmetry energy is almost 
beam-energy independent above 200 McV/nucleon. The 
same trend is observed for the single 7r~ /tt + ratio above 
1 GcV/nucleon. The sensitivity of the soft symmetry en- 
ergy to the ratio parameter is strongly affected by the 
choice of times, which is not true for the stiff symmetry 
energy. In simple words, just like the n~ ratio, the 
neutron-to-proton double ratio from free nucleons can act 
as a useful probe to constrain the high-density behavior 
of symmetry energy. Experiments are planned at MSU, 
GSI, RIKEN, and FRIB to determine the high-density 
behavior of symmetry energy by using the neutron-to- 
proton ratio. 
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